(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 43258, 1146] NotebookOptionsPosition[ 39979, 1025] NotebookOutlinePosition[ 40323, 1040] CellTagsIndexPosition[ 40280, 1037] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"c1", " ", "=", " ", RowBox[{"{", RowBox[{ "0.560", ",", " ", "0.560", ",", " ", "0.496", ",", " ", "0.496", ",", " ", "0.670", ",", " ", "0.670", ",", " ", "0.590", ",", " ", "0.590", ",", " ", "1.100", ",", " ", "1.100", ",", " ", "0.940", ",", " ", "0.940", ",", " ", "0.760", ",", " ", "0.760", ",", " ", "0.880", ",", " ", "0.880", ",", " ", "0.668", ",", " ", "0.668", ",", " ", "0.793", ",", " ", "0", ",", "793"}], "}"}]}]], "Input", CellChangeTimes->{{3.512566464470544*^9, 3.512566540592898*^9}, { 3.51256657279574*^9, 3.512566700692055*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "0.56`", ",", "0.56`", ",", "0.496`", ",", "0.496`", ",", "0.67`", ",", "0.67`", ",", "0.59`", ",", "0.59`", ",", "1.1`", ",", "1.1`", ",", "0.94`", ",", "0.94`", ",", "0.76`", ",", "0.76`", ",", "0.88`", ",", "0.88`", ",", "0.668`", ",", "0.668`", ",", "0.793`", ",", "0", ",", "793"}], "}"}]], "Output", CellChangeTimes->{3.5125667021681395`*^9, 3.512570441319007*^9, 3.513434147921034*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"c2", " ", "=", " ", RowBox[{"{", RowBox[{ "1.635", ",", " ", "0.922", ",", " ", "1.273", ",", " ", "0.718", ",", " ", "1.864", ",", " ", "0.655", ",", " ", "1.200", ",", " ", "0.245", ",", " ", "1.864", ",", " ", "0.569", ",", " ", "2.259", ",", " ", "0.0726", ",", " ", "2.640", ",", " ", "0.357", ",", " ", "3.337", ",", " ", "0.404", ",", " ", "1.742", ",", " ", "0.056", ",", " ", "2.980", ",", " ", "0.681"}], "}"}]}]}]], "Input", CellChangeTimes->{{3.512566708711514*^9, 3.512566882066429*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "1.635`", ",", "0.922`", ",", "1.273`", ",", "0.718`", ",", "1.864`", ",", "0.655`", ",", "1.2`", ",", "0.245`", ",", "1.864`", ",", "0.569`", ",", "2.259`", ",", "0.0726`", ",", "2.64`", ",", "0.357`", ",", "3.337`", ",", "0.404`", ",", "1.742`", ",", "0.056`", ",", "2.98`", ",", "0.681`"}], "}"}]], "Output", CellChangeTimes->{3.5125668894268503`*^9, 3.512570443627139*^9, 3.5134341479700365`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"c3", " ", "=", " ", RowBox[{"{", RowBox[{ "1.468", ",", " ", "1.468", ",", " ", "1.142", ",", " ", "1.142", ",", " ", "2.570", ",", " ", "2.570", ",", " ", "2.071", ",", " ", "2.071", ",", " ", "7.827", ",", " ", "7.827", ",", " ", "5.975", ",", " ", "5.975", ",", " ", "5.617", ",", " ", "5.617", ",", " ", "7.028", ",", " ", "7.028", ",", " ", "4.496", ",", " ", "4.496", ",", " ", "6.340", ",", " ", "6.340"}], "}"}]}]], "Input", CellChangeTimes->{{3.5125669376236067`*^9, 3.5125670808588*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "1.468`", ",", "1.468`", ",", "1.142`", ",", "1.142`", ",", "2.57`", ",", "2.57`", ",", "2.071`", ",", "2.071`", ",", "7.827`", ",", "7.827`", ",", "5.975`", ",", "5.975`", ",", "5.617`", ",", "5.617`", ",", "7.028`", ",", "7.028`", ",", "4.496`", ",", "4.496`", ",", "6.34`", ",", "6.34`"}], "}"}]], "Output", CellChangeTimes->{3.5125670926814756`*^9, 3.512570446022276*^9, 3.513434147992038*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"c4", " ", "=", " ", RowBox[{"{", RowBox[{ "0.765", ",", " ", "0.765", ",", " ", "0.868", ",", " ", "0.868", ",", " ", "1.053", ",", " ", "1.053", ",", " ", "0.944", ",", " ", "0.944", ",", " ", "1.122", ",", " ", "1.122", ",", " ", "0.962", ",", " ", "0.962", ",", " ", "0.864", ",", " ", "0.864", ",", " ", "0.987", ",", " ", "0.987", ",", " ", "0.671", ",", " ", "0.671", ",", " ", "0.875", ",", " ", "0.875"}], "}"}]}]], "Input", CellChangeTimes->{{3.5125670962786818`*^9, 3.512567219919753*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "0.765`", ",", "0.765`", ",", "0.868`", ",", "0.868`", ",", "1.053`", ",", "1.053`", ",", "0.944`", ",", "0.944`", ",", "1.122`", ",", "1.122`", ",", "0.962`", ",", "0.962`", ",", "0.864`", ",", "0.864`", ",", "0.987`", ",", "0.987`", ",", "0.671`", ",", "0.671`", ",", "0.875`", ",", "0.875`"}], "}"}]], "Output", CellChangeTimes->{3.512567220787803*^9, 3.51257044800939*^9, 3.513434148013039*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"c5", " ", "=", " ", RowBox[{"{", RowBox[{ "0.160", ",", " ", "0.160", ",", " ", "0.168", ",", " ", "0.168", ",", " ", "0.195", ",", " ", "0.195", ",", " ", "0.107", ",", " ", "0.107", ",", " ", "0.314", ",", " ", "0.314", ",", " ", "0.269", ",", " ", "0.269", ",", " ", "0.226", ",", " ", "0.226", ",", " ", "0.251", ",", " ", "0.251", ",", " ", "0.191", ",", " ", "0.191", ",", " ", "0.226", ",", " ", "0.226"}], "}"}]}]], "Input", CellChangeTimes->{{3.5125672428550653`*^9, 3.512567389146433*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "0.16`", ",", "0.16`", ",", "0.168`", ",", "0.168`", ",", "0.195`", ",", "0.195`", ",", "0.107`", ",", "0.107`", ",", "0.314`", ",", "0.314`", ",", "0.269`", ",", "0.269`", ",", "0.226`", ",", "0.226`", ",", "0.251`", ",", "0.251`", ",", "0.191`", ",", "0.191`", ",", "0.226`", ",", "0.226`"}], "}"}]], "Output", CellChangeTimes->{3.51256740342925*^9, 3.512570449963501*^9, 3.5134341480340405`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"s", "=", "s"}]], "Input", CellChangeTimes->{{3.5134348033469696`*^9, 3.5134348089161797`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"rezim", " ", "=", " ", "16"}]], "Input", CellChangeTimes->{{3.512569206881401*^9, 3.5125692154108887`*^9}, 3.5134341450288687`*^9}], Cell[BoxData["16"], "Output", CellChangeTimes->{3.512569217220992*^9, 3.512570452864667*^9, 3.5134341491731057`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"OmegaA", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"c1", "[", RowBox[{"[", "rezim", "]"}], "]"}], "*", RowBox[{"c4", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], "+", RowBox[{"c2", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.51256922557347*^9, 3.512569226423519*^9}}], Cell[BoxData["1.1280780114867943`"], "Output", CellChangeTimes->{3.512569228773653*^9, 3.51257045553682*^9, 3.5134341491951065`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"TA", "=", RowBox[{"1", "/", "OmegaA"}]}]], "Input", CellChangeTimes->{{3.512569237180134*^9, 3.5125692381201878`*^9}}], Cell[BoxData["0.886463515658825`"], "Output", CellChangeTimes->{3.5125692423714314`*^9, 3.512570458630997*^9, 3.513434149219108*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Tv", "=", RowBox[{"1", "/", RowBox[{"c4", "[", RowBox[{"[", "rezim", "]"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.512569251914977*^9, 3.5125692528540306`*^9}, { 3.5125705220676255`*^9, 3.512570545158946*^9}}], Cell[BoxData["1.0131712259371835`"], "Output", CellChangeTimes->{ 3.512569253438064*^9, 3.512570461143141*^9, {3.512570529014023*^9, 3.512570561713893*^9}, 3.5134341492451096`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Kc", "=", RowBox[{ RowBox[{"c3", "[", RowBox[{"[", "rezim", "]"}], "]"}], "*", RowBox[{ RowBox[{"c4", "[", RowBox[{"[", "rezim", "]"}], "]"}], "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"c1", "[", RowBox[{"[", "rezim", "]"}], "]"}], "*", RowBox[{"c4", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], "+", RowBox[{"c2", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], ")"}]}]}]}]], "Input", CellChangeTimes->{{3.5125692617505393`*^9, 3.5125692626595917`*^9}}], Cell[BoxData["5.450930407996479`"], "Output", CellChangeTimes->{3.5125692635306416`*^9, 3.5125704643433237`*^9, 3.513434149270111*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ash", "=", RowBox[{"2", "*", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"c1", "[", RowBox[{"[", "rezim", "]"}], "]"}], "+", RowBox[{"c5", "[", RowBox[{"[", "rezim", "]"}], "]"}], "-", RowBox[{"c4", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], ")"}], "/", RowBox[{"c3", "[", RowBox[{"[", "rezim", "]"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.5125692653907475`*^9, 3.5125692791435347`*^9}}], Cell[BoxData["0.040978941377347755`"], "Output", CellChangeTimes->{3.5125692796895657`*^9, 3.512570467211488*^9, 3.5134341492951126`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"bsh", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"c1", "[", RowBox[{"[", "rezim", "]"}], "]"}], "+", RowBox[{"c4", "[", RowBox[{"[", "rezim", "]"}], "]"}], "+", RowBox[{"c5", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], ")"}], "^", "2"}], "-", RowBox[{"4", "*", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"c1", "[", RowBox[{"[", "rezim", "]"}], "]"}], "*", RowBox[{"c4", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], "+", RowBox[{"c2", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], ")"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"c3", "[", RowBox[{"[", "rezim", "]"}], "]"}], ")"}], "^", "2"}], ")"}]}]}]], "Input", CellChangeTimes->{{3.5125692815586724`*^9, 3.5125692930413294`*^9}}], Cell[BoxData[ RowBox[{"-", "0.012234904596590464`"}]], "Output", CellChangeTimes->{3.5125692950574446`*^9, 3.5125704695506215`*^9, 3.5134341493191137`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mu", "=", RowBox[{ RowBox[{ RowBox[{"-", "ash"}], "/", "2"}], "+", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"ash", "^", "2"}], ")"}], "/", "4"}], "-", "bsh"}], "]"}], " "}]}]], "Input", CellChangeTimes->{{3.5125692965815315`*^9, 3.5125693075711603`*^9}}], Cell[BoxData["0.09200374246524261`"], "Output", CellChangeTimes->{3.512569308565217*^9, 3.512570471950759*^9, 3.5134341493421154`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Kcsh", "=", RowBox[{"Kc", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"mu", "*", "Kc"}]}], ")"}]}]}]], "Input", CellChangeTimes->{{3.512569312181424*^9, 3.5125693256351933`*^9}}], Cell[BoxData["3.630308781477998`"], "Output", CellChangeTimes->{3.512569326158223*^9, 3.5125704741618853`*^9, 3.513434149369117*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"TAsh", "=", RowBox[{"TA", "/", RowBox[{"Sqrt", "[", RowBox[{"1", "+", RowBox[{"mu", "*", "Kc"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.512569327542302*^9, 3.512569335596763*^9}}], Cell[BoxData["0.7234313588323567`"], "Output", CellChangeTimes->{3.5125693362367997`*^9, 3.5125704761069965`*^9, 3.513434149393118*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"If", "[", RowBox[{ RowBox[{"Tv", "<", RowBox[{"(", RowBox[{"10", "*", "TAsh"}], ")"}]}], ",", RowBox[{"{", RowBox[{ RowBox[{"iv", "=", RowBox[{"1", "/", RowBox[{"(", RowBox[{"Kcsh", "*", "Tv"}], ")"}]}]}], ",", " ", RowBox[{"nu", "=", RowBox[{"0.09", "/", RowBox[{"(", RowBox[{"Kcsh", "*", RowBox[{"Tv", "^", "2"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"iv", "=", RowBox[{"0.5", "/", RowBox[{"(", RowBox[{"Kcsh", "*", "Tv"}], ")"}]}]}], ",", " ", RowBox[{"nu", "=", RowBox[{"5", "/", RowBox[{"(", RowBox[{"Kcsh", "*", "Tv"}], " ", ")"}]}]}]}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.5125693376828823`*^9, 3.512569349507559*^9}, { 3.5125694013425236`*^9, 3.5125694064428153`*^9}, {3.512570415649539*^9, 3.5125704162745743`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.2718776995047141`", ",", "0.024150896047003746`"}], "}"}]], "Output", CellChangeTimes->{3.512569352174711*^9, 3.512569408605939*^9, 3.5125704802072315`*^9, 3.5125705366634607`*^9, 3.5125705680642567`*^9, 3.5134341494161196`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Tnu", "=", RowBox[{"iv", "/", "nu"}]}]], "Input", CellChangeTimes->{{3.5125694915476832`*^9, 3.5125695057204933`*^9}}], Cell[BoxData["11.257458065968706`"], "Output", CellChangeTimes->{3.512569515523054*^9, 3.512570484458474*^9, 3.5125705397616377`*^9, 3.5125705703433867`*^9, 3.5134341494401207`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"M1", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"s", "+", RowBox[{"c1", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"c5", "[", RowBox[{"[", "rezim", "]"}], "]"}], "*", "s"}], "+", RowBox[{"c2", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], ",", "0", ",", RowBox[{"c3", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"s", "+", RowBox[{"c4", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0", ",", "s", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"mu", ",", "0", ",", RowBox[{"nu", "*", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"Tnu", "*", "s"}], "+", "1"}], ")"}], "/", "s"}]}], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]}]], "Input", CellChangeTimes->{{3.512570632149922*^9, 3.5125706333159885`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"0.88`", "\[InvisibleSpace]", "+", "s"}], ",", RowBox[{"0.404`", "\[InvisibleSpace]", "+", RowBox[{"0.251`", " ", "s"}]}], ",", "0", ",", "7.028`"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"0.987`", "\[InvisibleSpace]", "+", "s"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0", ",", "s", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0.09200374246524261`", ",", "0", ",", FractionBox[ RowBox[{"0.024150896047003746`", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"11.257458065968706`", " ", "s"}]}], ")"}]}], "s"], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.512570635214097*^9, 3.5134341494661226`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"M2", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"s", "+", RowBox[{"c1", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"c5", "[", RowBox[{"[", "rezim", "]"}], "]"}], "*", "s"}], "+", RowBox[{"c2", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], ",", "0", ",", RowBox[{"c3", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"s", "+", RowBox[{"c4", "[", RowBox[{"[", "rezim", "]"}], "]"}]}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"mu", ",", "0", ",", RowBox[{"nu", "*", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"Tnu", "*", "s"}], "+", "1"}], ")"}], "/", "s"}]}], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]}]], "Input", CellChangeTimes->{{3.512570638847305*^9, 3.512570639202325*^9}, 3.5125706974896593`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"0.88`", "\[InvisibleSpace]", "+", "s"}], ",", RowBox[{"0.404`", "\[InvisibleSpace]", "+", RowBox[{"0.251`", " ", "s"}]}], ",", "0", ",", "7.028`"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"0.987`", "\[InvisibleSpace]", "+", "s"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0.09200374246524261`", ",", "0", ",", FractionBox[ RowBox[{"0.024150896047003746`", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"11.257458065968706`", " ", "s"}]}], ")"}]}], "s"], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5125706992597604`*^9, 3.5134341494921236`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"W", "=", RowBox[{"Simplify", "[", RowBox[{ RowBox[{"Det", "[", "M2", "]"}], "/", RowBox[{"Det", "[", "M1", "]"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.512570701491888*^9, 3.5125707125515203`*^9}, { 3.513434247017702*^9, 3.513434263265631*^9}, {3.5134343049130135`*^9, 3.513434306042078*^9}}], Cell[BoxData[ FractionBox[ RowBox[{"1.9107564721191301`", " ", RowBox[{"(", RowBox[{"0.08882999999999996`", "\[InvisibleSpace]", "+", "s"}], ")"}], " ", RowBox[{"(", RowBox[{"0.9870000000000003`", "\[InvisibleSpace]", "+", "s"}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{"0.09819131525378877`", "\[InvisibleSpace]", "+", "s"}], ")"}], " ", RowBox[{"(", RowBox[{"0.8439203931775128`", "\[InvisibleSpace]", "+", "s"}], ")"}], " ", RowBox[{"(", RowBox[{"2.021657564112605`", "\[InvisibleSpace]", "+", RowBox[{"1.8224905936144233`", " ", "s"}], "+", SuperscriptBox["s", "2"]}], ")"}]}]]], "Output", CellChangeTimes->{3.5125707314095993`*^9, 3.5134341495201254`*^9, 3.513434309941301*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"s", "=", RowBox[{"I", "*", "om"}]}]], "Input", CellChangeTimes->{{3.51343416324191*^9, 3.5134341636639347`*^9}, { 3.5134342895171328`*^9, 3.513434314054536*^9}, {3.5134346705907364`*^9, 3.513434737639654*^9}}], Cell[BoxData[ RowBox[{"\[ImaginaryI]", " ", "om"}]], "Output", CellChangeTimes->{{3.5134346783751507`*^9, 3.513434743754865*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData["W"], "Input", CellChangeTimes->{{3.5134347464380703`*^9, 3.513434747670472*^9}}], Cell[BoxData[ FractionBox[ RowBox[{"1.9107564721191301`", " ", RowBox[{"(", RowBox[{"0.08882999999999996`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}], " ", RowBox[{"(", RowBox[{"0.9870000000000003`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{"0.09819131525378877`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}], " ", RowBox[{"(", RowBox[{"0.8439203931775128`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}], " ", RowBox[{"(", RowBox[{"2.021657564112605`", "\[InvisibleSpace]", "+", RowBox[{"1.8224905936144233`", " ", "\[ImaginaryI]", " ", "om"}], "-", SuperscriptBox["om", "2"]}], ")"}]}]]], "Output", CellChangeTimes->{3.513434749433275*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"F1", "=", RowBox[{"Abs", "[", "W", "]"}]}]], "Input", CellChangeTimes->{{3.5134347564688873`*^9, 3.513434771366914*^9}}], Cell[BoxData[ RowBox[{"1.9107564721191301`", " ", RowBox[{"Abs", "[", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"0.08882999999999996`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}], " ", RowBox[{"(", RowBox[{"0.9870000000000003`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{"0.09819131525378877`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}], " ", RowBox[{"(", RowBox[{"0.8439203931775128`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}], " ", RowBox[{"(", RowBox[{"2.021657564112605`", "\[InvisibleSpace]", "+", RowBox[{"1.8224905936144233`", " ", "\[ImaginaryI]", " ", "om"}], "-", SuperscriptBox["om", "2"]}], ")"}]}]], "]"}]}]], "Output", CellChangeTimes->{3.5134347729113164`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{"F1", ",", RowBox[{"{", RowBox[{"om", ",", "0", ",", "20"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "1.2"}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.513434820663*^9, 3.513434877993101*^9}, { 3.513434930284393*^9, 3.513434962202049*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwVlHc41e8bx41IZH8JycjKjpMicj9ICA0rSUgiGeEgq5IRQllFygyVmUJ2 x8jeI7tsdTicT0hWfuf313O9ruu57ue+7/f7eQvb3jW6RUNFRYWjpqL6/5ny wJZewe6V+l1hClCRobFfyOCKkC9MsxDkCdRk2GuqDRYSigDBmIE0Ij0ZWOLK PjAOx8HJTFMysJIhhRAtSCf0Ek542Xy04iODzIrd0z2HNLhixzfdKUGGCZbU fztDmbBLWHvKcIoMzme8XDYE34Pn9JOCHnMyfBqYzFrTyodjj07rut0hw5aT wTjmUAR0ozsGBwLIcLYzVuzP0EcIHbUtYEwlw88knbo1xxLYuLRQwFpIhii7 XcvV7VJw0k7U0qohQ//27QSyYAXEu5y5sjtMBu9mgeMrxZVwVrAg1n+ODHzx A+0krWpwtWyZ5sTIUGv1xGFpsAZsYqI3OrfJYCuFaBcdvoCgNKHpEx0GNNt+ 3WsBBFAhYu96WDAYmDQ9sqZdBwF6ebnHD2Hgk89U8XuoHm6LVSwTxTA4HzdP /zu9AdSXJQ9cl8GA36fOBHNshIr2W081FDFYvv46k6z4FXKdYlnenMKAoHWP vLL9FUpVLJ0j1TCIkzRSX2lsgscyoiO0GhjYscpGLUc3A3Nw77f/tDE4ub5/ lGTWAu7vQs806WKwf2xagiTYCoHB5VqCBhgME2q8ln62gsr9ujDxixjk5iQ1 LBa3wc0iTr7pyxgEROHZF/3a4al4ftgVEwwueFywJmp1QJSKxWykGQZC5pIF vw52Ai/3u4Fgcwx+n9m3/XOwEx54Q6KWBQaNIj90f6Z2gQFnGXP3NQxeHKh8 seDQDRL25Y2S1zEouXfeA8/TA52kmG0zKwxeE0/xrgb0AI0Ga95VawxCr4sR PKZ6oCB6MPKEDQauPRwOv7V7wUKUi59I4StaVCweub1A39qQFXADA1RGKsFY +sB83Vh9mcKSkmPX3PF9QK22NYBsMeB43UKDDfXBfGmEjDeFt1nK3rup9YPQ InnmGYVnH725RE7vB7NDNj7/5861mI27dANgzCCs+//7ZQ4PUlccB2BD/thz bQqnjTpp3+0agOrvt8L+Ud4LN7y6tKw4CCa7E3KZFHYnnIt3TRwEY/qoUAUK W+BOnF7eHgRtTx73Akr/WjnCUy4234BkZtfITWEZXtZwUuM38NBfQa6U+bmi duRcJIdAM+dwUgllP//+/Rpcih4CJVUTZyJlfwvuQwHOv4eg7Bu3ERuFe2Yb RZbMhiFEar5H3BKDjLY090XBEVgNnPKXoOixllQq9SNqBBYZfZ5tUfQ659A+ 07c5AuJOGeE1VzBYot0wrewfhV/HNqv+M8UA+g6yFqIxyPJoSn9vTPFb+tGW jIIxqFxoqZM3wkD5zIXTEWHjEKR6OJmN4p9IJrvVgLVxGHgznetkiMH3Ed98 txsTMCzSnFeuj0GId7aAuep3eDfDNnCM4seeoh0aiZUfkEWW8rdFFH3ggBWz 3CTYrvPsu6COgap4wukwu0mIX3jmIU3x+9O196t+vZMg0iMYUkX5DydiB2/Z 5k0Bf88IJi2HgZWPjeb41BRwz/3RyZTGIMx6UcCMZxryTPkZD0hiMCpLPawX Og2ipHlCjggGD9tkzitYz0DQf6fz+XgweFf8WTzv+Qzc+DKDO8iFQV+SJq1Y xwyMfJB6t8iOgaiDeTWPyiyIXE0Y82DCoI02RI6Kcw4MjuEfJ+6SYY3IwuSn NwfTafLppE0yHOl7ubD6cA6uSs96yv0hg1t6UfrC0hy4rBTGuS6TgfvMGEd3 0zxMnbbCjU+QAYnar+juzsOpoRuN/iNkcGTC2utxC5BSfiKAZpAM1SP0oWXp CyA2//ZJQzslf7wVN1L8fsIXhw8X95eTobAoYsxZjgi3ymDheCQZiBzEqPM3 iGDbw7ffO5QMEt7n1Y8lEIH5Lqt92kMypKsxZcxsEuF1sHNGCp4Mca1R9hZf F6Gly6IBLMjQJUM6pPx3Eexy1UJTjcnAFGPYyiW9BJNpC/DDgAwhZiwyvTFL cBb9y6MDSn7OPMN0LEkw0HPbYlKYDB/PkTPFnpGAqK7ZrkvJ+5X3l0xo60lQ 9n7ydSwHGW67sX+ulViGGe2FgQYaMljsxgUorS6DXdgWw8WGFUhjag9m16b0 2bqe/UVxBdTP+FsGWlH4APd/kYIrMO4qrbRyjwyJ1crG3AdXgKcvcr4jlww1 3huiLrPLEJtkoBfOhkFxYBCLaNwyBIt1seyNYfBlSGmCfZoEDtD7csljFVoO Vl6yMF8Cqsw/lWORq5D1KNk0R30Jkvbxj7dlrcKZOg3DVtElaGm1F3j/bRXY C1x5ksiLcMxkJ/PW6TUYx10heT1ehJ+3JfK/U6/DzT+eDsL5RLgTd7+2J+YP 8O5IvCXN/wTXOcmZksJNWFojyIWZzMGk1QhDXesmlA29drE8OQfGw+FynbOb cFVZjk2PZw6U2xd85/i2oD9P2dppfBZoi7PZuMO2QHdOi37h1iwkBQir37Pe hnYH5x54MAN1nLwvVdh2IZW/+Wty8xSciG6p1ZbeBafK9iirginIofeZvXxu F9q7XY/pxE9B1OY3+TsBuyDiIOIeZj0FZj8SmpJ/7sKlc+KvmzcngZjLtrZN +AetN6sVrihNAqNwiMQPZSpk/+XFZ+fYCWhRi5V/oE2FCowrmeUdJ+Cxeeqp I0ZUKFqvVVRBYwJoYsp1rjlRIeGMwIZVbBy2/i05DKVQIU0H8Sxns3FYHDd9 101Djewf+NBtSI1BR5KEJKGDGhXe3DwZRDUCT0pOKFiPUCMfw1/5NT+GQbdH Q+XfHDViWl09ePzLMDTut9RT26NGRfcrNcsfDEO1d4xjuQIN+uA35W9INQwF JpvvP7ygQWNK7raazEPwlK1dKsOGFvXuMBayGgxCxakjpRdcaZGUUoPcZflB mLW6Czv+tGgrZPtAM8cgqBZwmlxJpEU5DpanL4wOwE89y/vMXbSIpP015IvT AGgFk7p9Vfch+05526HEfvi7zuplxE2H3M6wqGlz9IHIEdu9PRE61EB3U+DK di9cOFsSUaBAh0rcEgtiZ3shO848lcGQDnlQ9wU//NwLRvIZzbXBdOijiiCj knUvFNxW5JPG6FB5T4JRXWkP2I6ZfKHtoEffr78/nrzVBYyvplzcR+jRDIeR psO3Lvho4cr/Y54eKUewtbp97ALa0TDfSur9yNFv4vwppy7IGq7EuZ/ajyIQ thjyoxPmB4Xefc/cj452RepV93XAnZ6lmApfBmQ6ZdKQPNkGHDG+IBHGgBQz nrAdbmiDyov0ywkJDEhSedKwLbsNmLoF9d2KGFD1WI5IiVMb5Hca0UnMMSAe jVFc3mYrLLeV+yZcOoB+1J1RHTnSCu5fQ2zvSjCi5jBSynJIMwjMnqP/g2NE u5d9ZXVcm6Gd5kBuAGJEaXaHi75caQYxFI1FXGVELNV33rBKN8NI5fPArEhG 9BLFi4n1N4HGh5zUkRVGtFc6WfLnWBOwv2oZ1SpnQqnDQ8/afzVCbcWTB22N TKhfpC+j6lsjOA0bHL3cy4RiE4QNvzY0wleuXkcrIhNq7hisFU1tBN+Ykb8+ /AeRkWmfPqtJI0yFErkLAw8iKZui/1wbG6DY7aAxjx4zOvtsN121rB46kjl9 C8yY0f7kSeFPWfUw38iXpmnHjN567f7WjK8HPl7JReeHzKiJf99Qrls9hNRp B9eVMiOZp6m/X8jUgxnHw093jrKgvdf0dNXv6mDrE8ZZs82CbP8YPZcoJsB/ 3/+qGB1gRcfiOJcsMgkgx0Bls8DNij7uShAy4glw05Iln12RFfkQL1s4ehOg Y5/0WfvbrGhBZrtqXI0AqaZ2nqyDrMhbinTPLfkLZPg7EY9fYUPjDjxz39dr 4DjhwN+w62yI++PmdfPvNUDY947ux00KE3UsSE01MBk9KxTtxoZa8m1OOCbV gED69Su/ItjQVRvJ6TdqNZDceLExo5oNdZrLVlM/qYY45hOpHEfZkdzDdBbj U1Vw1Kgvz/EYOwpqUP1IPloFxS/cKghy7GhDqftINksV9AgWDLiqsqNDsZNN l+crgUVRnKndhB0pWU82lL2ohCdmPD7BYeyIyZZJK367AoLSdi6vLbGjdCnS U8JAOVw6Xu4uvcqOBG+ZDxg3lINAHT7WdpMd8WcrWG8Wl0PlNLGnh44DtZdw T/o9K4ffYiMXCwQ40FnlWwMm+uVgm19qaH+JA1Gf2TlzsekzaFS46g1/4kCB qY82qlrKgOW8lCNrJQe6xl207lhRBuOjc+HnCBxop63MUiK3DHx2LFtKOjiQ M40dZ09kGXwAfZ3YOQ60+Blb+naxDISaJLTPH+JEJHrswPvRUqDqn0TVfpwo YZmHb2yrBEJ+1Wb4PeJE1/CqfVJLJXCAKoVGJZwT7c6cXImYKAEO2auNpc85 Ua3VZRpvQgmIhfXqFH3gRPSzakAbVgLnVesvZM5zokn1Cd5J7hKIf/PGMtzo P2TBd9aTT+MTiHva+5hIcaH4WI3ZoqpikD39wzjtOBeaevUZz1BcDCeozOWJ J7nQ5fETr+7kFINmlN78Qy0uJObvFKIfWwxWWTImeZZcSJcjMuiBQzG8GPgt T/OMC80vUPmlcRUDndLDhaJVLlTwVXw7yOcDzK4nmjHVcKMoQ/dhi4JCyNTA 4ScauJGO7lHq5uRCsI7uelbUxo0WIjCGs+GFMCpK12o8zI0EhQLt3G4WQq+x h+rrVW40dJOPJ4avEGqLDYRkpQ4heZZu5dHIAkh0oSZeSDyE7kYFeHH55cP5 uTv349x5kNx6XKTVo1w4d18k/6wPD9Kgszxb4ZELmlzjo38e8KBnzlbCQna5 cFrbUPlaFA+SeBDCwaGTC1LZ8qsib3mQ7bUoq1DmXGC6tXa7dJwHvS0S10x9 /R46Zu+bDOvwors0fvMaDe/gwmyMtIAAH+o37Alrk38LHdKzA5aifEji252B FNG3oI8/9eCVFB8y+SB50pf3LejSfO/hOcWHSvrXeHVp34KmkLQ35yU+RKog PdEeygEly68EhiA+lPR20L08MAf4BzZNV+f4kHT/lID+SDYQG2wDWwoOI9sn XXey07JAMHb4ztOSw0ih/4eb//MsMLG6YGpSdRiV+8Zsm0dmQe1fFanJlsNo jwv3SOpeFsTJsg9szBxGddPygeIXs0A18YuEBB8/Ypnr3fOnyoKoO/zdoY/5 kVm0vrXz7Tcgx/7tiJb1EdT3e6RqQz8TGGH9xviqABqzq/jwwjsddEgWG0dN hRB3gkI3TWQKDNfbUVvICCO6526cd0uSoTS80UpoXRhdPmFlJOGYBPZ0XpNd XUfRu2W9fj71F3DTw7/0YLIICmvtkef1SQB5B+kugxuiKDCDW6qDFAfdHDJM XipiCPc09/VNj1ioKgx6JUMrjoYWgtalT8WAP7XUpY4JcaSxb2wym+kpMI77 6bqWSCDGqasEl9VImMsRTkgLOIaCpVzeXjz+BBKO8e6rMJVEQXayDmVe4RCT bN98Q0gKaUz8Nd/69xgsEq53HVyQQl+y9Eysi0Lh79nszbRaaST5ssS+2DYE ZuY0ferDZNChtbDCmHPBkHlG4XDnDVnUMYuF0MgEQUGoF7uDrBwyZq31DVd5 BGGmN/Vp1+UQvkUurvF6ICiNZ7wP65RH87/T4zCah5Ak+aqgU/84ynKrV7J6 fB8eEqKT9CuOI88IdUVWhgAwPGya3SaqgIpcjvw5m+8HGq7N19gjFFA5nbeL +ANfSOO3Fe1YV0BiW3uyRQY+INrw9sWFa4rI9KvbJ2/NexAjWy/VV6uIthLo mqoMvcF+dV9vqxQOvfrUeHLBxAtsi38LRTzDoTDnw+58/p7gHGS8cz4WhwyE 0z473fMEL+OSoYPxOFSU5LBTg/eE8HWvpzEvcEhW8erDW86eUKiyuZ2YgkNb NS02hOuesFm39y0nD4fELPI1s5EnxPQdjG5sxqHDylbskvs9IfmNy+3HrTiU WmXl8ZLWE954dmnptuPQZndwEyOVJ5Rxx2y1d+FQPDZ1jvwXD+MWnLf7B3FI zzEotYmIB4kZXq3pGRyqrcBbpXfh4XiJn0DWHA7deBreKtiOh9OhY5u3FnBo TN/5aHozHgwkUj78IuLQG7VLr98Q8ODhJCSAYThUr2DztvQjHgLUHm1+XKXU D6m7p1aEh1Dm6QHPdRwK+WuCa8zDQ1LRm8i/fyn1UprufMvCQ2bgPofKLRxS Ntubs8nAQ97lW5oBO5R51AX1l1LwUHK06Yj6PxwqXe9/6ZOMh9pV8c29PRy6 pt7eS5eIh/8BNmQ9VQ== "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, ImageSize->{958., Automatic}, PlotRange->{{0, 20}, {0, 1.2}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Automatic}]], "Output", CellChangeTimes->{{3.5134348516758547`*^9, 3.5134348788199024`*^9}, { 3.513434933685199*^9, 3.51343496291965*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Integrate", "[", RowBox[{"F1", ",", RowBox[{"{", RowBox[{"om", ",", "0", ",", "10"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.513435045662195*^9, 3.513435142819166*^9}, { 3.513435251410957*^9, 3.5134352546245623`*^9}}], Cell[BoxData[ RowBox[{ SubsuperscriptBox["\[Integral]", "0", "10"], RowBox[{"1.9107564721191301`", " ", RowBox[{"Abs", "[", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"0.08882999999999996`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}], " ", RowBox[{"(", RowBox[{"0.9870000000000003`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{"0.09819131525378877`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}], " ", RowBox[{"(", RowBox[{"0.8439203931775128`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "om"}]}], ")"}], " ", RowBox[{"(", RowBox[{"2.021657564112605`", "\[InvisibleSpace]", "+", RowBox[{"1.8224905936144233`", " ", "\[ImaginaryI]", " ", "om"}], "-", SuperscriptBox["om", "2"]}], ")"}]}]], "]"}], RowBox[{"\[DifferentialD]", "om"}]}]}]], "Output", CellChangeTimes->{3.5134351524287825`*^9, 3.5134352584465694`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"h", "=", RowBox[{"InverseLaplaceTransform", "[", RowBox[{ RowBox[{"W", "/", "s"}], ",", "s", ",", "t"}], "]"}]}]], "Input", CellChangeTimes->{{3.5125707338827405`*^9, 3.512570734288764*^9}, { 3.5125708062898817`*^9, 3.512570808415004*^9}}], Cell[BoxData[ RowBox[{"1.0000000000000002`", "\[InvisibleSpace]", "-", RowBox[{"0.27430444514612395`", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "0.8439203931775128`"}], " ", "t"}]]}], "+", RowBox[{"0.11721255984221624`", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "0.09819131525378877`"}], " ", "t"}]]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "0.9112452968072117`"}], "-", RowBox[{"1.0914621263055084`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "0.42145405734804625`"}], "-", RowBox[{"0.45263957911433916`", " ", "\[ImaginaryI]"}]}], ")"}], "-", RowBox[{ RowBox[{"(", RowBox[{"0.42145405734804625`", "\[InvisibleSpace]", "-", RowBox[{"0.45263957911433916`", " ", "\[ImaginaryI]"}]}], ")"}], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"2.182924252611017`", " ", "\[ImaginaryI]", " ", "t"}]]}]}], ")"}]}]}]], "Output", CellChangeTimes->{3.5125708092920537`*^9, 3.5134341495581274`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{"h", ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "1.2"}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.512570810543125*^9, 3.512570829199192*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV13c4lt8bAPDXeIeVPZPRq5QiRb4NuW+KJGlSSstIRaIdSoNIKQlpIirJ yMxKr5G9V0JRZI/3Edn1O7+/nutzPc91znnuc8597qNqe3qXAy+NRrvFQ6P9 /8nDXVG/YnEo0DzHV0tz9cHOiXlFw/Qx0OLtbgXWGkFKQ0f02MZncK3G79Rn j80w7WTeRjmGA+00V5tOmUNizrXRWIFI4EytDAnh3QnrZtirl9KiAKX8nlLz LaE3bHPe2IloQB4TNnf3PrhrP2fze+YVXPtvfd0uDRvQ1k6ZpALeQMTKpp59 7CPQ3BUu/EXoLdAaHW0u/GcLzuPvYWFLLHQcXe1j/8seFB42lA9tjIMIGfOq 8QWOYKuBfAOOidDxs9HXx9UZeGfcq8c830OEuGtAgvlpaOiwXDBmnAScr0uP GS5yg0txQpmjX5LBtfiVq57lOTAL6maMRqQAN5B7TpxxARQv5e2hTqTCjlkP U9vii8DZeJE7MpMGYh0VMgwPD2C2/lQfUs6EwBvBwu2116GZ8/H8YG8m1Gw1 9PnSfwNiX4cVDCRlQYzUX4wT8gaLMxaH+zfmgGn+xOts+1sQKpAV2uP4CcQS RcOWnb0DqRfNzpyV48ASr/nHP5XehWf9/8n/9uQARzTNuU79HrjUSDiOGudB tEYSn9qVQJB4VsJLfcmHigh/izCPhzAzL/2tq34BHFm63stgUzB0XY/awY0o ABWPxIcXxEIg3fHqi5EThdBLccp9UkJhv47uuuGZz6Bb2DVyg/0EIsvC3QaU S0Gbtf6MUH0EjIWlabTfLQVPZYFBp4ORYOJY3lk3VQrz6F4j7/siYZBvwjKr vgzkLmZmtDCjYM0Gi3W3fStgIa/OD/P9r6AmcZZXfaQadN/uHk7cEwteZcvN Vh6uB4eVeysEC5MhJunD4nch9bCC/boR1FKgLsyIb1FFPTg92iexyicF1Bz3 5citbQDDkC9fNm9JhTI+by2aZCMsn78/Q7glDWQ2tEpUFzWB72t/zUjFTEhI vN3qrNUCDK3DpUkiHOiX6L9rdrQFGr5n/DXaxAH1C2YGS4JbYNVerSpxDw5E 6AtFdk61gMJ05JBDHweCSu8e2/+5FX7fG1vaXJIHFzrvU5ttvoEFtVrqWFAB hAuV3xQ37oDEO7IJ28yLwRFqHw+e6YTFObnOCwyqweXX0s7UhD5QMv+ovHfq Cwiqequ3r+HCsvPdKx+EdsA9sXKNyCO/oXDJ249Fb7rAtnXPJ76KP9Bc4fF3 6HEvCD79ccrt6x8wPZNQbRLTC8n7XRTbu/+A4fMK3dD0XuBr8b2cxTMBqjND b+TqeyG6OUvH7b8JSB+dSHon1AfdjSox319OQN/75BCGRx+crBkMzLw8CWEa P8+c39UPbp+9bU+rT8PWw/rD4Z0DoNRlwvijMw06m4o4H4YHoJxXINYTp8EL 4uuKpgZgEQZQt62nISmkN7lMdBC+ZoVci74zDWWMogMn1g+C4fvXL76OTENN hqo674NBEH9a0rIxYwbWZm11f7l6CHIz/a+WFc7AX7e9mi9hCJyazRfurJ2B kvKvix5vGYLP0rUnDvXPwM92v37ng0NwOfDr5CXFWXDTPtv6znsIfvj0yyRc m4WNc0fAr3YIklyFd8ttmQO5dXsHCuyHoeKJ5OV4qzl4WF363sVlGLoLFcKN 7OdAe7LUXvTSMCjILx1w9poDdbf2yFV3hsE7z/hmXtoc7FoSkEpLGgYrCa+U kwv/wsNf3658nB6G6RRK8uPMX9C/IBtReWsEpL5Prt0l8A8kDtpfXX9/BLRY tCM9MsT5yTHhj0bAzmZenPiqf+Ces3XELGYEKviXbTp2/B/89Rx5o1k6Ai8s 7c+JNhIfX3z1IT8XIj2c+rX30vCUdqwmryEXtDkCk74HaeiwSETYfhMXOPwx 9HY7Gk48jOnO28yFjoAulQBXGr72inzkZsEFpYiDe/tu0zBOqkgh8gAXnhRu L4zMoWHIpLWB0XkuBInovpBYyIN34os0nrziwsJdde9OLOHBUy46S7xiuJAU 6prJ0eJBwXsBh46840KNcnyDy3oe/JileFU2iQvzVi0WKt/Dg6ZGPnn7crjg byV36aYvD6rdClq0tI4LN8Jnd44N8mB1tpWo2AwXdmhnuC37zYMvEs5TZ+fI +PLOPrCdIu1L/TjT8I8LWT/7a2rovLh+rvXAHX4KRhd93R6vxIs/s3bQakUo sI1L23ZsBy+GH852+KhMgWGmy5bmFF685bLGbqkhBfPMNE6IZvHiqyJFg0Ub KWhr+eVnwuHFGoPOmQXGFFyatSlJreBFK8eMJYwtFLyHrZsf/OLF9rp53kk7 KFApUjc2k+VDn4CT6hcOU0Cr78Acdz4sWibtWOROgXdfbqT7dT489bc138eT AgHac961fnxoT8soMrxKgYSmdWFaCB9qO1G1idcpWORbuznxPR/Oym40tvej wGx9vsXLbj40uZ4hszKEgqqdEYlHhvjQeW2fcU4oBbuPXxVTHuPDtkZTs01h FNiErKt7ysOPOzK/fTN7SoHLSLJliCI/pu044LUqkoKHUVE2frv48Y+ZOdsl jgK5rOsfTaz5cZzZOFEcT8GzmsNK9CP8mGr/ukopkYLXc/M7rp/ix7Eou3sF SRRk7A228/DjxwXL/0b3pZN4CfmcdMnlx8A6u6hYDgWLzx27tEeDjp98tlyQ rKNAc1377nBtOvb3Ws/y1FOgS9u3ol+Pjm0de1YPExvd3dLttZGOLCta6adG Cg5FL9/zzoaOXsplaZu/UhDaMLqC9z4dxdVswgbbKXj+xEloWwgdo8Pq933o oCD6SFf3o6fEsjUjXj8oSB5sfK4ZQ0f+hctiBTtJvPgzhazz6PgrSugeo5sC +mqvnsTfdFwqeG7u8QAFwjMT+dNTdFw+8Vpz4yCZjzzXF8Y0Bt6LO7K3n1hl m61lqzADT+wUj1g1TIG+g3EBYzEDh259M3rHpeBciFD4wX0MpG93TDQYp8Dj gLd7zCEGJjhvcy8hvqE6Z/nbnoH6dNbRHX8oCIwfFr7txsB26bP++ycoiPtc 657qz8B99gZ++6Yo6Bp/ZCX0kYEDYR30uVkKXhrqnP1WwMDPtteDzsxRcDig 6n5iGQOX8NeKdBO3qNFLdzczUEiYp7XoLwW1u8+sf/abgRnds3KutFG4Fy6y z2WagVefn7dvJDYfiDmHPEw8fdInZQ3PKJTcaI/vmsfECY+iY1PEuUnmKpoa TNzxS+ikHd8oeM716P/TZuLIw5sKmcRrt9y0rv2PiWpKSzuF+UchtSMz6Lwx E7kG2sUJxLGi6vTcI0zsX2qs2kYfheMH8lUDHZkYeHxry0LGKCx6c9DA1oWJ 4Jv3yZE4wiD4IsOTiUEP+RUGiR+d4um3eMTEr5yyne3MUdiT+Yyh+oKJRicM UmVZoyBBX8P+Hc1Ew1XCey2I7z1zsXmUzMTgZ4GO6cTmvQKXT2Qy0XRAu6WH WED3Vch6DhM5yX9DZQVGwbu8taq9kokRm7sn3YiNZC8OJDUw0dr0+8tnxP9s JVjerUxkLeZEfCZ2nzY1XNLHxLK+oUgJwVFYY9J1cHqEiQtthWP1iMcfeLlX /GFi4tdMUWti16XpKa78LLyQbzUXRqx1fmeNkRALu/4J+6cTD3AGB6UkWJhT 136ljjhG2E+gR46FB/q1WgeJj+1jL85UJt9f3BdNFxoFtehcozuLWainmdmq SPxzxPrwQU3yfn3czVXEB30Dw3jWs1AsOkHLmlihfllavSEL2U3OeieJm5WK a1+ZsrBEbH7GZeLQk7bDF7ezUGlcON2XeE/6nKCZFQtlXwSvCiaW4H2srniQ hX7hvNrhxNXbdDcN27FwA/31+xjiu4+rj3BOslByV0X6e2KzXyevBLmxUM4r b9MHYpMr7LhNl1goobTdNIfYSLqt5c9VFqoIx/36RGwQHyzw1oeF7w1oyvnE 64y3rTlwl4Ww7syf/1vvG91R5CELw8TVXf7vVedzQz49ZqG8g5Efh1hL5GKh WwQLp492m3wk1ni14jf7DQtD8nTiMogXb+hVbYpnYUrJlqxk4oWNETv8UlnY ZmJy9h2x0ilrr3XZLIwPMf768v/xokskDOaxcEmY7XAYsczzsrYXJWT8Wemp Af+Px+qbQjurWbh1xFrrOvG8yvXr+JpYuDLF9uBZYiGHseNpbSycY3ZvtCdm zcU9cuxkoeaFic7dxPwhDkXy/SzsJncMI+LZwib21QkWOnhprppPPGlzf5f2 XxY+FltQ+f/5HhvbfP0nvwBKr+heNkLWw6Ba5ncTCQFkKGltyCXuzXETmZIT wNaMj7+jiLv2aOi/UxbA1HOmF/yI27yfPhbVFMBE1yM924gruq7sad4sgAsO 5CaVkfVc4ql3099CALOtLGlviAulRpL0LQWwPfqnxg3inE1HRCNtBfCTVt78 1cRx0UZlJ64IYOzQ4agAsp9i9Gcm53sL4CK+W6yjxNENKepV/gIYwNW31iF+ xr/IZ1WYABZniHxuIPs3wJ6JM8kC6Hw8Y0aQ2H/2k0t8pgBeNh8abST7/Vbw peeHOQL4eslgVzjx1cK+6YJKAbxpUczRJnZRq0i72yuAu1UV7pqQfGLRFbhM SUkQG0cPiWwi+ahiWVeDjZog/tlC6fATbz3739WnGoIYhCmBBbyjYMr7vUbu P0Fc2nxnFIiNVJZdkNwhiAZj6V4rSL5bbfOZw7ohiJ5f98/1k/yp2DBl+fuX IC7/m6/zjeTjJ/O3/Vs5KIi2a0WabhHL20XEuI4Kop2kZ8sKYplRk5mhv4J4 SWNy3rVJCsTEgiN6ZYVw3VVPM3mS3/nMtQa+bRHC7AceS7THKOgvsL1WEi+E egEFSyeHKFB+0HzyXqoQipg5p/gS7zlkYbknWwi97XOTZIhzJ9dqdJQIYfJT lNAh51GQpnjDRKcQnlkTdMehn4L1jz6pqysI46GpgYhkcr7dPalY7XNLGKdX vBCb+E6BlnjTgo2HRbC0+1yZWxUFgjB+tO33PFTI+7TYn5z3m4f2Tyy0FMOa VgjJvk1Bc749z/7l4pjzRyqg+zgFaX6Fh1TGxdFY8cmJxG0UxJiVy/lOiqNX 8tN8Q+InwnV1QzPi2PC1z6HBnIJrge0m2bwSuOx63Z3prRSYh01rWYlJ4JKX 1/eamVHw8432vzvLJFBKtJBn0oTUX8XPwieOSqBawA+lu0jBMfr5jqoqCfR6 ePjFY10KFJRLv8zUSqC6OG1mE3HVmgXVSxolUK5t0xOuDgV6zp8/3mgl3ytF R5gSM+qln+r1SiDQ6RUzKyl4FZFu+YJHEqWK34a6rCDn6/rJslO6kvjzY2y6 jwYFdmc80oSfSGKW6U2pi6qkvVXf05Y+l8RD88SkdYkFRzHdJEISv9Ro7qJU SP3nxvhw7bUk1omXXHcmnnV9kDGWLIlCd5VTHUh9GXr6TXZbuSSuY/HtsltA QalzfV7cnCRuyG+v8pcn9dVyvfxSmhR2Hz+1ejux62BYfjefFM7bFtooSSzj fKhARVAKXz/Ykf9cjtSvTn2FwTJS2HX7qlmqLOnvxL9izxVSaMER+DsgTcEK x2VV5kelkGbtoukjQcH3E8qn4+2lkHM5doE58V1nSbF5x6XwAo/vNkniPreZ ndUuUlgyfsvppTgFUZ7ljTs9pTDIw/ZZgRgFskFO36weSWHrLrlEcVEK/uW8 GzxSKYWRvIFXW4UoSPgUHpBXI4UtfO+lY4lt8h9qLWyQwoK8oz2XiDOLPVw7 W6Rwt1/tWlnis3Vbxxx6pVBy0iDUSpCCnp6BmZN80thycpTZwaKgWmK50Pm1 0jizZDRckUGBwalTfU/1pXF2R6rAOJ30V5xQnA/SqHc5P7KS+J7nSm8xE2lU zV5l40Vs3q03926XNLrXuP7pIvV/WabhyE8naSzVzYzI4KPg85G99TteSOPb psHKazyk/swOS7oYKY02b3uzDhFHS7fcfxEtjQ+somv1ib3LbMwHY6WxKbvE Y4pG6lNd28++H6TRqbgg9xwxh3nqQ26NNMZsdms/Q+4f2Qk3ni7nk8G4+tV3 7pL7CcPcUuERQwadWiqmnIl39qk/5hGUwZlvgsHbiHvYlaFNYjIosyHfTJRY 8rFs0DUlGTxySDgneJYLTjfj/BrWymDeRLDPG3L/mb+v6bynqwzOddidGZji wrHxmLHuszI4bHRjtoY4Kcjj7M6LMsg55R+XTry5UsVt8VUZ9NsQdeUG8Vkj J+eaOzKo6cVsn09cvpxmp/ZGBn9u3nbRcpILHjwaOyq+yaBLzamin3+40L1J qHnxDxnMmORTqiTe4Td4+HqXDC6fqgz+QKwmmnhab0AG+4qt2+8SVyjq3o+Y lEHRTS2da4iV1hhUnZOUxYBEum/oOBfyXXaZL9gii69f8yo5jnFhebJOw0Vz WfSzkp9nSRw6LmVTt10WtSeOy20kdvL84uRrJYvzKuU8lIml/G3uUHayuOud 77mvv8n/vzpW9vmqLL5PtJDZRSzY5m7qkiqL/jn7RS1HuaB38a3I/Q+yGLV3 O20z8VGJ5rrELFl8HsbiX0ucYbr6IMWRRcXmPRsXEB9LG3E7VymLYgUaL7op LuTds3/q3i2LamlZ7zyJLxlaDN+Sk8MRRtjRNC4Xots8U97Ml8OO+uXn3hJX X3x3qURJDo30maHPiRcnsPgEF8nhl2ZV5i3ieoUC2YCVcii8lv5sL7HW2H+G D83kUL/C/fLcCBd+vVYNDveUwxdn8y12Ew96/5jiuyaHKzZsX29KPGobefj4 TTlcclNIbwMxTVl12Up/OfxwNNRSnVjhkUp+/iM5DHAIXTs7TObLV5n7K1kO TR9EuMcQ5zguMF/eJ4cvQz3FBYmDl8jzZ1rKY2XujOvAIBdsDi1+07lPHl82 6Yd0EqsF65jNs5FH/8K6T63Eqf+2BdrZyqPJ1gT1CuKGphuK807LY/0zGsYT S/oM6tr5yePKpetuuBIH/fjkIJItjzphTmXTA1wIfHKs+KiKAmoeF3mkQLw/ +GCVcI8CymVeo5/v48LkpldT4bnz0WUqW+dJDxc6fxldyvdVxMPGA7pxv7jw csPK+ZVHF6DD49wPJp1ciPc5L+6oqYRrbSf0Mjq44Gtpt5VvXAmrL7y80f+N C6vbIt/6ViqjZ1ewQlALF8KWPo2v3KqCKz7+09P9wgUvTkDY1kwV5DTzS+2r 58K2+ZavytRU0fjWPd5j1VwwdCk+IH5bFVNHzcojyrkQrmirVjGuigu6GIFX ikl8Ct6EWhxYiNY9/GEvCsj/aeZr1OUuxNsOg2+bP5H19Ju/tlSDjSxjb4fK bC7YJo2q3L7Pxk8H+QtkM7jgfGP3rNkDNmq3Jh9gEp/fnfpF+CEb4y+2jP/5 wAW/8fP3AkPZSKmEaTURJ6ydmnn0nI079/9MCiaeyvvX9PodG7/xWNEkiAPr hAMKi9loGeaVJ5LOhSdRp47fKmVjrHui/1waF6LOVW00LWdjVO8pyyHidJnA 6fIqNua9DRupIG7bL3m8vpGNaevS9O4Sq3fKb/zZycbHZ66LCRFrp7orRf8i /Se2jE6ncmGdT+uUQw8bq2vLmvqJzdWfv+/rZ+NsrFtUGfEZJxUlimJj+4Pu 7f7EnvrXp5J/szGpnw/ciX1EfjacG2djV2n2ipPEYYlRdyYn2dhRMyttRvzy Gr9j1jQbe6pvCa0jfrfTwchzlo01fPG8GsSpC4sWGPxl4+eG0zPyxLm/F0/9 +8fGh2+rxwWI/wfsMZ8d "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0, 10}, {0, 1.2}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Automatic}]], "Output", CellChangeTimes->{3.5125708308782883`*^9, 3.5134341495931296`*^9}] }, Open ]] }, WindowSize->{1350, 633}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (November 10, 2008)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 613, 11, 31, "Input"], Cell[1183, 35, 454, 9, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[1674, 49, 602, 11, 52, "Input"], Cell[2279, 62, 461, 9, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2777, 76, 557, 10, 31, "Input"], Cell[3337, 88, 458, 9, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3832, 102, 559, 10, 31, "Input"], Cell[4394, 114, 459, 9, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4890, 128, 559, 10, 31, "Input"], Cell[5452, 140, 459, 9, 30, "Output"] }, Open ]], Cell[5926, 152, 119, 2, 31, "Input"], Cell[CellGroupData[{ Cell[6070, 158, 160, 3, 31, "Input"], Cell[6233, 163, 120, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6390, 170, 386, 11, 31, "Input"], Cell[6779, 183, 136, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6952, 190, 146, 3, 31, "Input"], Cell[7101, 195, 136, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7274, 202, 251, 6, 31, "Input"], Cell[7528, 210, 187, 3, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7752, 218, 564, 17, 31, "Input"], Cell[8319, 237, 138, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8494, 244, 485, 14, 31, "Input"], Cell[8982, 260, 141, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9160, 267, 957, 30, 31, "Input"], Cell[10120, 299, 160, 3, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10317, 307, 340, 11, 31, "Input"], Cell[10660, 320, 138, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10835, 327, 221, 6, 31, "Input"], Cell[11059, 335, 136, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11232, 342, 227, 6, 31, "Input"], Cell[11462, 350, 139, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11638, 357, 951, 30, 31, "Input"], Cell[12592, 389, 280, 6, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12909, 400, 146, 3, 31, "Input"], Cell[13058, 405, 185, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13280, 412, 1133, 35, 31, "Input"], Cell[14416, 449, 890, 24, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15343, 478, 1159, 36, 31, "Input"], Cell[16505, 516, 892, 24, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17434, 545, 342, 8, 31, "Input"], Cell[17779, 555, 759, 20, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18575, 580, 242, 5, 31, "Input"], Cell[18820, 587, 131, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18988, 594, 95, 1, 31, "Input"], Cell[19086, 597, 876, 20, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19999, 622, 148, 3, 31, "Input"], Cell[20150, 627, 958, 22, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[21145, 654, 355, 9, 31, "Input"], Cell[21503, 665, 6898, 119, 628, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[28438, 789, 264, 6, 31, "Input"], Cell[28705, 797, 1103, 25, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[29845, 827, 277, 6, 31, "Input"], Cell[30125, 835, 1211, 31, 38, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31373, 871, 307, 8, 31, "Input"], Cell[31683, 881, 8280, 141, 253, "Output"] }, Open ]] } ] *) (* End of internal cache information *)